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ABSTRACT 

Polarizability tensor of a strongly magnetized plasma and the polarization vectors and opacities 
of normal electromagnetic waves are studied for the conditions typical of neutron star atmospheres, 
taking account of partial ionization effects. Vacuum polarization is also included using a new set of 
fitting formulae that are accurate for wide range of field strengths. The full account of the coupling 
of the quantum mechanical structure of the atoms to their center-of-mass motion across the magnetic 
field is shown to be crucial for the correct evaluation of the polarization properties and opacities of 
the plasma. The self-consistent treatment of the polarizability and absorption coefficients proves to 
be necessary if the ionization degree of the plasma is low, which can occur in the atmospheres of 
cool or ultramagnetized neutron stars. Atmosphere models and spectra are presented to illustrate the 
importance of such self-consistent treatment. 

Subject headings: magnetic fields — plasmas — stars: atmospheres — stars: neutron — X-rays: stars 



1. INTRODUCTION 

In recent years, thermal or thermal-like radiation 
has been detected from several classes of isolated neu- 
tron stars (NSs): radio pulsars with typical magnetic 
fields B ~ 10 12 -10 13 G, "dim" NSs whose magnetic 
fields are mostly unknown, anomalous X-ray pulsars and 
soft gamma-r ay repeaters with B possibly ~ 10 14 -10 15 
G (see, e.g.. iBecker fe Aschenbachl |2003 lHaberll [2004 ; 



llsrael. Mereghetti. fc Stellall2002HPavlov fc Zavlrtl2*003 . 

for reviews). The spectrum of thermal radiation is 
formed in a thin atmospheric layer (with scale height 
~ 0.1-10 cm and density p ~ 10 _2 -10 3 g cm -3 ) that cov- 
ers the stellar surface. Therefore, a proper interpretation 
of the observations of NS surface emission requires un- 
derstanding of radiative properties of these magnetized 
a tmospheres. 

Shib^ngv^t^y (^22) (see also iShibanov fc Zaviinl 
1995: Pav lov et al.lll995l) presented the first model of the 
NS atmospheres with strong magnetic fields, assuming 
full ionization. Variants of t his mo d el were con s tructe d 
by iZane. Turolla. fc Treves! (|2000ft; iZane et alJ (|2001l) : 
IHo fc Lail J200l 12001 : lO^el ll200ll) : ILlovdl ll200l . An 
inaccurate treatment of the absorption due to free- 
free transiti ons in strong mag netic fields in the ear- 
lier models dPav lov et a fTl 19951) has been corrected by 
iPotekhin fc Chabrierl l)2003|k this correction has been 
taken into a ccount in later m odels l |Ho et al Jl2003l 120041: 
I Llovdl l200l . Recent work l|Ho fc Lail l2003t lLai fc Hoi 
120021 l2003a|) has shown that in the magnetar field regime 
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(B > 10 14 G) vacuum polarization significantly affects 
the emergent spectrum from the atmosphere; for weaker 
fields, vacuum polarization can still le ave an unique im - 
print on the X-ray polarization signals l)Lai fc Hol2 003b'l . 

Because the strong magnetic field significantly in- 
creases the binding ene rgies of a toms, molecules, and 
other bound states (see lLail 120011 for a review), these 
bound states may have abundances appreciable enough 
to contribute to the opacity in the atmosphere. For 
calculation of this contribution, the non-trivial cou- 
pling of the center-of-mass ( CM) motion of the atom 
to its internal structure (e.g.. lPotekhinlll994l and ref- 
erences therein) can be important. Also, because of 
the relatively high atmosphere density, a proper treat- 
ment should take account of the plasma nonideality 
that leads to Stark broadening and pressure ioniza- 
tion. Recently, thermodynamically consistent equation 
of state (EOS) and opacities have been obtained for a 
magnetized, part ially ionized H plasma for 8 x 10 11 G 
< B < 10 15 G iPotekhin Chabrier. fc Shibanovl H999T 
IPotekhin fc Chabrierl2003l l2004D . These EOS and op ac- 
ities have been implemented bv lHo et alJ 1)20031 120041) for 
modeling NS atmospheres. For the typical field strengths 
B = 10 12 -10 13 G this modeling showed that, although 
the spectral features due to neutral atoms lie at extreme 
UV and very soft X-ray energy bands and therefore are 
difficult to observe, the continuum flux is also different 
from the fully ionized case, especially at lower energies, 
which can affect fitting of the o bserved spectra . For the 
superstrong field B > 10 14 G. IHo et alfl|2003ft showed 
that the vacuum polarization effect not only suppresses 
the proton cyclotron line, but also suppresses spectral 
features due to bound s pecies. 

It is well known (e.g.. lGinzburg|ll97Ci IMeszarosll 1992ft 
that under typical conditions (e.g., far from the reso- 
nances) radiation propagates in a magnetized plasma in 
the form of two so-called extraordinary and ordinary nor- 
mal modes. The polarization vectors of these modes, e x 
and e° are determined by the Hermitian part (x H ) of 
the complex polarizability tensor (x) °f the plasma. Our 
previous treatment of these modes in partially ionized 



atmo s phere s i|Potekhin fc Chabrieill200ll2004HHo et all 
120031. 12004ft was not quite self-consistent, because the 
effect of the presence of the bound states on the polar- 
ization of normal modes was neglected: we adopted the 
same pola rization vectors as in the fu lly ionized plasma, 
assuming (|Potekhin fc Chabrierl f2003') that the effect of 
bound states on these vectors should be small provided 
the ionization degree of the plasma is high. However, this 
hypothesis (related to % H ) may be called into question, 
based on the observation that the absorption coefficients 
(corresponding to the anti-Hermitian part of the complex 
polarizability tensor, ix A ) are strikingly affected by the 
presence of even a few percent of the atoms. 

In this paper, we study the polarizability tensor, the 
polarization vectors of the normal waves, and the opac- 
ities of the partially ionized nonideal hydrogen plasma 
in strong magnetic fields in a self-consis tent manner, us- 
ing th e technique applied previously bv lBulik fc Pavlovl 
( 1996) to the case of a monatomic ideal hydrogen gas. 
In fJ2]we introduce basic definitions and formulae to be 
used for calculation of the plasma polarizability (new fit- 
ting formulae for the vacuum polarizability are given in 
the Appendix). An approximate model based on a per- 
turbation theory, which explains the importance of the 
CM coupling for plasma polarizability, is described in SJ3 
The results of numerical calculations of the plasma polar- 
izability are presented in 21 an d the consequences for the 
polarization and opacities of the normal modes are dis- 
cussed in JO In fJHl we present examples of NS thermal 
spectra, calculated using the new opacities, compared 
with the earlier results. In ^TJwe summarize our results, 
outline the range of their applicability, and discuss un- 
solved problems. 

2. GENERAL FORMULAE FOR POLARIZATION IN A 
MAGNETIZED PLASMA 

2.1. Complex Polarizability Tensor 

The propagation of electromagnetic waves in a medium 
is described by the wave equation that is obtained from 
the Maxwell equations involving the tensors of electric 
permittivity e, magnetic permeabilit y u, and electr ical 
conductivity <r. It is convenient fe.g.. lGinzburdll970|l to 
introduce the complex dielectric tensor e' = e + Airier /uj. 
In the strong magnetic field, it should include the vacuum 
polarization. When the vacuum polarization is small, it 
can be linearly added to the plasma polarization. Then 
we can write 



e' = I + 4ttx + Attx 1 



(1) 



where I is the unit tensor, x = X H + *X A is the com- 
plex polarizability tensor of the plasma, and x vac is the 
polarizability tensor of the vacuum. 

In the Cartesian coordinate system x'y'z' with unit 
vectors x', y', z', where z' is along magnetic field B, 
the electric permittivity tensor of a plasma in the dipolc 
approximation, the dielectric vacuum correction, and the 
inverse magnetic permea bility of the vacu um can be writ- 
ten, respectively, as fe.g. JHo fc Laill2003l and references 
therein) 



I + 47TX 



e igO 
-ig e 
77 

diag(a, a,a + q), 



(2) 
(3) 



fi 1 = I + diag(a, a, a + m), (4) 

where a, q, and m are vacuum polarization coefficients 
which vanish at B = 0. The formulae for calculation of 
these coefficients are given in the Appendix. Quantities 
e, 77, and q are well kn own for fully ionized ideal plasmas 
fe.g.. lCxinzbu73ll970Th In this paper we shall calculate 
X for partially ionized hydrogen plasmas at typical NS 
atmosphere conditions. 

The complex polarizability tensor of a plasma becomes 
diagonal, x = diag(x+i, X-ij Xo) in the cyclic or rotat- 
ing coordinates, where the cyclic unit vectors are defined 
as eo = z', e±i = (x' ± iy')/^/2. The real parts of the 
components \a ( a — iljO) determine the Hermitian 
tensor x H , which describes the refraction and polariza- 
tion of waves, and their imaginary parts determine the 
anti-Hermitian tensor ix A , responsible for the absorp- 
tion. According to equation 

4vrx±i = £ - 1 ± g, 47TXo = »7-l- (5) 
Note that in the cyclic representation, the general sym- 
metry relations for the polarizability tensor take the form 

X A hu) = - X A a (u), x2(-w)=X?». (6) 

2.2. Relation Between the Plasma Polarizability and 
Absorption 

General expressions in the dipole approximation for \ A 
and Xa through frequencies and oscillator strengths of 
quantu m transitions in a magn etized medium are given, 
e.g., bv lBulik fc Pavlovl l)1996ft . For transitions between 
two stationary quantum states i and / with energies 
Ei and Ef — Ei + hwfi > Ei and number densities 
of the occupied states ni and nf, the absorption co- 
efficient for the basic pola rization a equals (see, e.g., 
lArmstrong fc Nicholislli972ft 

27r 2 e 2 

tfl = pk a {uj) = (rii - n f ) Pj5{u - u fi ), (7) 

m e c 

where f*f = (2m e u; fi /h)\{f\r\i) 
tor strength for the transition i 
tial absorption coefficients sum up into the total fi a — 
Si f(E f >Ei) tfa j where Jw/ includes integration over 
continuum states. Then the equation 



1 2 is the oscilla- 
/. These par- 



Xa{uj) = 4^ qH ' 



(8) 



together with the first symmetry relation JSJ) yield 



— y 

i,f(E f >E i ) 



fj5{u-u H ) + fl a 5{u + u fi )\. (9) 

Once x A is known, x H can be calculated from 
the Kramers-Kronig rela tion llLandau fc Lifshit3 119891 
§123) or its modification ijBulik fc Pavlovll 1996ft 

1 poo 1 A/ f\ 

x^)--v ^M^d.', (10) 

7TW J_, 



where V J means the principal value of the integral. From 
equations (JSJ and (jTU)l . 

e 2 fif fif 

He \ e 1 \ J a J — a 



2 m e oj 



OJfi 



UJ fl_ 

(11) 
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Taking into account equations (JSJ) and © , we can present 
relation I jlOjl in the form convenient for calculation at 
uj > 0: 



x2H = 



4n 2 u 



duj' 



fi a (uj +U 1 ) - fJ, a (oJ ~ T 



jU a (o/) 



du/ 



duj'}. (12) 



If we replace delta-function in equation J7J by the 
Lorentz profile {v x J[ fir) [ (uj — Lo a ) 2 + (^J) 2 where 2/^/ 
is a damping rate for a given transition, then equations 
(jHJl and (|11J) can be combined to give 



2 m„o; ^— ' 

i,f(E s >E % 
if 



(m - n f ) 



J a 



if 



J —a 



OJfi 



if 



(13) 



3. ANALYTIC MODEL OF POLARIZATION OF ATOMIC 
GAS: EFFECT OF CENTER-OF-MASS MOTION 

As mentioned in §T\ in strong magnetic fields the in- 
ternal structure of an atom is strongly coupled to its 
CM motion perpendicular to the field. This coupling 
(referred to as "CM coupling") has a significant effect 
on the radiative opacities and dielectric property of the 
medium. Before presenting numerical results for the po- 
larizability tensor of a partially ionized plasma in 21 it 
is useful to consider an analytic model to illustrate this 
CM coupling effect. 

Consider an atomic gas in which the atom possesses 
only two energy levels, with the upper level having the 
radiative (Lorentz) width i/, assumed to be constant for 
simplicity. The energies of the (moving) atom in the 
ground state and the excited state are denoted by E\ (P) 
and E2(P') respectively, where P and P' are the CM 
pseudo-momenta, and the subscripts 1 and 2 specify the 
internal degree of freedom of the atom. [In strong mag- 
netic fields, the internal quantum numbers are (s, v z ), 
where s = 0, 1, 2, • • • measures the relative angular mo- 
mentum between the proton and electron, and v z is re- 
lated to the number of nodes in the z direction.] If 
there were no CM coupling, we would have E\^(P) = 
£1,2(0) + P 2 /(2m H ) and uj 2 \ = (E 2 - Ei)/h ^constant. 
Therefore, in this case equation (|13fl would yield 



XaM 



e 2 n 1 
1m e uj 



■I a 



f 12 

J —c 



U21 



■ W21 



( 14 ) 

This would imply that even for very small neutral atom 
fraction, the bound-bound transition can severely affect 
the dielectric property of the gas in the neighborhood of 
uj = uj 21 (for v -C W21). 

However, in strong magnetic fields, the energy associ- 
ated with the transverse CM motion of the atom cannot 
be separated from the internal energy. We therefore have 



e 2 m f Pl {P±)d 2 P ± 



2m e uj J (2irh) 2 
uj + uj 2 i{Pj_) + iu 



/i 2 (Fx) 



w 2 i (Pi 



(15) 



where P± is the transverse pseudo-momentum, 
Pi(P±)d 2 P± is the probability to find an atom in 
the initial state "1" in an element d 2 P± near P± 
(see iPotekhin et aO 11999^ . and we have used the fact 
that the oscillator strength is nonzero only when 
P = P' (in the dipole approximation). The CM 
coupling effect will significantly smooth out the di- 
vergent behavior at uj = uj 2 \ in equation (|14|) (for 
v — > 0). This effect can be taken into account usin g 
the perturbatio n result of [P avlov fc Meszarosl |l993), 
as was done bv lBulik fc Pavlovl (Il996f>, or n umeri cally, 
using the techniques of | Pavlov fc Potekhinl \ 19951) and 
IPotekhin fc Pavlovl l|1 9971) . 

In the perturbation theory, valid for small P±, the 
coupling reveals itself as an effective "transverse" mass 
AI± > mu acquired by the atom, which is different for 
different quantum levels. For the two-level atom, the 
perturbation theory gives Ej(P±) = Ej(0) + P|/(2M J _ i ) 
(j = 1,2), and we find 



Xc 



e 2 ni 
2m e uj 



dy 



fl 2 e~ 



•W2l(0) 



■a(T)y 



12 



■W2l(0) 



a(T)y 



(16) 



where a(T) — (kBT/h)(l — M± 1 /M±2), and we have used 
Pi(P±) oc exp[~P 2 /(2M ±1 k B T)}. At v ->• 0, only weak 
logarithmic divergence is present near uj = W21, com- 
pared to the much stronger (uj — W21) -1 divergence in 
equation l|14|l . Additional line broadening due to col- 
lisions (which implies v ^ 0), thermal effect and pres- 
sure effect, not treated in this simple model, will further 
smooth out the "divergent" feature in the polarizability 
tensor. 

This simple model shows that a proper treatment of 
the CM coupling is important in calculating the polariz- 
ability tensor of a partially ionized plasma. In practice, 
since the a bsorption coefficient iin has alread y been cal- 
culated bv IPotekhin fc Chabrierl pOOl 12004]) . we find it 
is more convenient to apply the Kramers-Kronig trans- 
formation (eq. ^21) to obtain Xa ( see ©• 

4. EFFECT OF PARTIAL IONIZATION ON 
POLARIZABILITY OF A STRONGLY MAGNETIZED 
HYDROGEN PLASMA 

Some approximations for calculation of the complex 
dielectric te nsor of a p l asma have been discussed, for 
example, bv IGinzburd l|197fJ|) . In particular, the "el- 
ementary theory" gives the expression widely used in 
the past to describe polarization properties of the fully 
ionized electron-ion plasma (e.g . . | Shibanov et alJ Il992t 
IZaufi et a,l]l2?KKllHo fc Lall2n?fll l2f 



47TXa 



(17) 



(uj + auj ce )(u ~ auj c i) + iuju e g ' 

where uj ce = eB/m e c and uj C i = ZeB/m^c are 
the electron and ion cyclotron frequencies, uj p \ — 
(47rn e e 2 /r7i e ) 1 / 2 is the electron plasma frequency, v c s is 
the effective damping frequency, m\ is the ion mass, Ze 
is the ion charge, and n e is the electron number density. 
For the hydrogen plasma, the corresponding energies are 
huj ce = 11.577 B u keV, huj ci = 6.3049 B 12 eV, where 
B12 = B/10 12 G, and hu pl = 0.0287^/p keV, where p 
is in g cm~ 3 . In general, the damping frequency v c r 
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0.2 0.3 
hcj (keV) 

Fig. 1. — Absorption coefficient (upper panel) and polarizability 
(lower panel) of a hydrogen plasma for longitudinal polarization 
(a = 0) at B = 2.35 x 10 12 G, p = 0.1 g cm' 3 , and T = 3.16 x 
10 s K, according to four different models: fully ionized plasma 
(dot-dashed lines), partially ionized plasma without the effect of 
coupling between CM motion and internal structure of the atom 
(short dashes), with the magnetic broadening taken into account 
by perturbation theory (long dashes) and numerical calculation 
beyond the perturbation approximation (solid lines). 



depends on polarization and photon frequency. Deriva- 
tion of equation (|17|1 is based on the assumption that 
the electrons and ions lose their ordered velocities (im- 
parted by the electric field of the electromagnetic wave) 
at the rate v e s which does not depend on the velocity. 
A more rigorous kinetic theory gives Xa ano - Xa which 
cannot be in gene ral described by equation (11711 using 
the same v s (e.g.. lGinzburdll970l §6). For instance, if 
we describe Xa using equation l|17jl with some frequency- 
dependent u e fi(u)), then the Kramers-Kronig transforma- 
tion will give Xa which, in general, does not coincide 
with the real part of equation (|17|l with the same v e g(uj), 
although, for realistic i/ e g(u>), the difference should be 
small atw> w p i. 

Realistic models of the absorption coefficients of a par- 
tially ionized plasma do not allow us to perform the 
Kramers-Kronig transformation analytically. We evalu- 
ated the integrals in equation i|12|) numerically. Since 
the integrand can be sharply peaked near resonance 
frequencies, we e mployed the adapt ive-stepsize Runge- 
Kutta integration ijPress et al.l ll996. §16.2). The accu- 
racy of the numerical transformation was tested with the 
models that allow one to perform the Kramers-Kronig 
transformation analytically - those given by equations 
EH) and 113 above , and by equations (49)-(53) of 
Bulik & Pavk>3 <11996l) - and proved to be within 0.1%. 

Let us first consider a model that neglects the CM 
coupling and the plasma nonideality. The absorption co- 
efficient of the fully ionized component of the plasma in- 
cludes contributions due to the free-free absorption and 
scattering on free electrons an protons. For the atoms, 
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Fig. 2. — Same as in Figure HI but for a = +1. 



we include in consideration the bound-bound transition 
with the largest oscillator strength (for every polariza- 
tion) and the bound-free transitions. In this case the 
absorption by the atom can be described by analytic for- 
mulae. The bound-bound absorption cross section is de- 
scribed by the Lorentz profile, where the effective damp- 
ing width is mainly contributed by t he electron impact 
broadening ijPavlov fc Potekh in|| 1995 1) and can be evalu- 
ated using a fitting formula (Potekhin 1998). The bound- 
free cross sections are well described in the adiabatic 
approximatio n (however this is not tr ue with the CM 
coupling, see lPotekhin fc Paylovlll997D: thev are repro - 
duced by fitting formulae in iPotekhin fc Pavlovl |l993). 
For B = 2.35 x 10 12 G, p = 0.1 g cm" 3 , T = 10 5 5 K (the 
neutral fraction xn ~ 0.12), the resulting absorption co- 
efficients are shown by the short-dashed lines in Figures 
E] and for a — and +1, respectively. For a = this 
absorption profile is similar to the idealized model con- 
sidered above. The corresponding x^^) are shown by 
the short-dashed lines on the lower panels of Figures E] 
and El In each figure, dot-dashed lines correspond to the 
model of a fully ionized hydrogen plasma at the same B, 
p, and T. As in the above analytic model, we see that the 
presence of the atoms results in strong deviation from the 
fully ionized plasma model in the vicinity of the principal 
atomic resonances (the bound-free resonance for a = +1 
is almost invisible on the lower panel of Fig. El because 
of its small oscillator strength). 

The absorption coefficients obtained using the pertur- 
bation theory of magnetic broadening (Q are shown by 
the long-dashed lines on the upper panels of Figures El 
and El The long-dashed lines on the lower panels show 
the corresponding polarizabilities. We see that the reso- 
nant features are smoothed down by the magnetic broad- 
ening, and the resulting curves of Xa( w ) do n °t much 
differ from the fully ionized plasma model. 

Still greater smoothing occurs for the accurate (non- 
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Fig. 3. — Polarizabilities v? (a = ±1,0) of the partially ionized (solid lines) and fully ionized (dot-dashed lines) plasma with the same 
parameters as in Figures H|2| The right part of the figure has an enlarged vertical scale. The dotted lines correspond to uj p \ (vertical) and 
(q + m)/(2n) (horizontal), where q and m are the vacuum polarization parameters. 



perturbative) functions /i Q (w) and Xq( w ) (solid lines in 
the figures). The absorption coefficients shown on the 
up per panels correspon d to t he opacities in Figure 9 
of iPotekhin k. Chabrierl l)2003|) . The significant differ- 
ence of these absorption coefficients from the perturba- 
tional ones arises partly from the dependence of the os- 
cillator strengths on the atomic velocity across the field, 
partly from transitions which were dipole-forbidden for 
the atom at rest and were not taken into account in the 
perturbation approximation described above, and partly 
from the nonideal plasma effect of continuum lowering. 
Nevertheless, despite these differences in p, a , the differ- 
ence in is much less significant. 

Figure [3] shows all three components x^ m a wider 
energy range. In this figure, as well as in the previ- 
ous two, the polarizabilities in the fully ionized plasma 
model, shown by the dot-dashed lines, are obtained from 
equation 1|12|) . where the free- free contribution to /x a (w) 
includes the frequency-depen dent Coulomb logarithm 
ijPotekhin fc Chabrierl 12003'!. Because of this frequency 
dependence, these polarizabilities are not identical to the 
ones given by equation 117(1 . the difference being small 
at uj ^> cjpi and appreciable at lo < u> p \, where both 
the elementary theory and the description of absorption 
by binary collisions (i.e., neglecting the collective mo- 
tion effects) are rather inaccurate. The solid lines, which 
are obtained for the partially ionized plasma, are fairly 
close to the dot-dashed lines. The only prominent fea- 
tures are the proton and electron cyclotron resonances 
at huj ci = 0.0f48 keV and huj ce = 27.21 keV. 

However, such smoothing does not always occur. Let 
us consider a higher magnetic field B = 3 x 10 13 G 
and density p = 1 g cm -3 , retaining the same temper- 
ature. Then xh = 0.89. Figure 01 demonstrates the ab- 
sorption coefficients and polarizabilities for the fully ion- 
ized (dot-dashed lines) and partially ionized (solid lines) 
plasma. In addition to the proton cyclotron resonance 
at frui = 0.19 keV, absorption coefficients in the partially 
ionized plasma exhibit magnetically broadened bound- 



bound and bound-free features. The most prominent are 
the bound-bound absorption at fiw w 0.2-0.3 keV for 
and the photoionization edge at ftio — 0.408 keV for 
/io. These features are clearly reflected in the behavior 
of x^+i an d Xcn shown in the lower panels. Thus, with 
increasing B, the abundance of neutral states increases 
along with their influence on the polarizability. 

A similar trend occurs with lowering T. If, for exam- 
ple, in the previous case (B = 2.35 x 10 12 G, p = 0.1 
g cm~ 3 ) we take a lower T = 1.58 x 10 5 K, then we 
will have 94.1% of H atoms in the ground state, 1.3% in 
the excited states, and 1.1% of H2 molecules. At p = 1 
g cm~ 3 , the abundance of the atoms will be 96.9% (with 
roughly equal fractions of excited states and molecules), 
and w e nearly recover the case studied bv lBulik &: Pavlovl 
(1996), who assumed 100% of atoms for these plasma pa- 
rameters. 

5. EFFECT OF PARTIAL IONIZATION ON POLARIZATION 
AND OPACITIES OF NORMAL MODES 

In the coordinate system xyz with z along the wave 
vector of the photon and B in the x-z plane, the polar- 
ization vectors e° of the n ormal modes in a magnetized 
plasma can be written as (|Ho fc Laill200l 



(p3 f)3 jdJ 
4i °y 1 ^zj 



(i 



K^ + Kl^^K^l^K^), (18) 



where j = 1, 2 correspond to the extraordinary mode (X- 
mode) and ordinary mode (O-mode). The parameters 
Kj and K z j are expressed in terms of the components 
of the dielectric and magnetic tensors as 



Ki = p {1 + (-iy 
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where 



r)' - e' + g 2 /e' + rj' m/(l + a) e' sin 2 
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Fig. 4. — Absorption coefficients (top panel; the lines marked 
"+", "— ", and "0" correspond to the polarization index a = +1, 
— 1, and 0, respectively) and polarizabilities x+i (middle) and 
(bottom panel) of the partially ionized (solid lines) and fully ionized 
(dot-dashed lines) plasma for B = 3 X 10 13 G, 
T = 3.16 x 10 5 K. 




1 g cm , and 



9b is the angle between B and the z axis, and (see 
Eqs. PP-jHl) e' = e + a, rf = i] + a + q. In the usual 
case where the plasma and vacuum polarizabilities are 
small (|Xal ^ (4 71 ") -1 and \a\,q, \m\ <§; 1), the polariza- 
tion parameter (3 is given by 



X+i 



X H i 



(q + m)/(2ir) sin 2 



2(x?i-X H i) cos b 

The opacity in the mode j can be written as 



^| e ^( w ,c?B)| 2 ^H, 

a=-l 



(22) 



(23) 



where k a (a — —1,0,1) do not depend on Ob- For a 
partially ionized, strongly magnet ized hydrogen plasma 
K, n (co) have been obtained by IPotekhin fc Chabrierl 

(20o.i l2nol . 

The polarization vectors e x and e° for the polarizabil- 
ities presented in Figure prove to be indistinguishable 
from the results for the fully ionized plasma at general 9b 
values. They exhibit vacuum polarization resonances at 
huj ~ 10 keV related to intersections of x^-i with the com- 
bination of vacuum polarization coefficients (q + m)/27r 
that enters equation H22J) (the horizontal dotted line in 
Fig.|3 in this case |xc?Ij Ix+il ^ IX-il) an d the electron 
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Fig. 5. — Squared moduli of the cyclic components e 3 a of po- 
larization vectors e x and e° for the normal waves propagating at 
angles Bg = 60° (two upper panels) and 10° (two lower panels) 
in a hydrogen plasma with B = 3 X 10 13 G, T = 3.16 X 10 5 K, 
p = 1 g cm -3 . Solid lines: partially ionized plasma; dot-dashed 
lines: full ionization. 



cyclotron resonance at tkj as 27 keV. If we neglected the 
CM coupling effect, we would observe additional resonant 
features near ftui — 0.07 and 0.2 keV, associated with the 
bound-bound transitions for a = +1 (Fig. |2J and a — 
(Fig.^). Actually these features are smoothed away by 
the CM coupling. 

For the second set of plasma parameters considered in 
21 the ionization degree is relatively low (1 — x-r = 0.11). 
Figure |S] shows squared moduli of two components of the 
polarization vectors, e 3 and e 3 +1 (the third component 
can be found as \&-\\ 2 = 1 — |e+i| 2 — | e 1 2 ) for two prop- 
agation angles 6b = 60° and Ob = 10°. Dot-dashed and 
solid curves correspond to the fully and partially ionized 
plasma models, respectively. At 9b — 60° (two upper 
panels), there are two sharp resonances for the partially 
ionized case at hoj = 0.158 and 0.425 keV, associated 
with the zero level crossings by Xoi^) (the bottom panel 
of Fig. 0}, which are absent in the case of full ioniza- 
tion. The feature at 3.3 keV is the vacuum resonance. 
At smaller angle (two lower panels of Fig. O, these res- 
onances become broader, and there appear additionally 
other features associated with the behavior of x+i (see 
the middle panel of Fig.0J). For the fully ionized plasma, 
an additional feature is just the proton cyclotron reso- 
nance, whereas for the partially ionized case the behavior 
of the polarization vectors is more complicated because 
of the influence of the atomic resonances. 

Figure shows the opacities calculated according to 
equation l|23() . using \e J a \ 2 shown in Figure [S] and k a (uj) 
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Fig. 6. — Opacities in two modes Kj (Eq. 1231 ^ for the same plasma parameters as in Fig. 0and two angles, 9b = W° (left panel) and 
8b = 60° (right panel). Solid lines: a self-consistent calculation for a partially ionized hydrogen plasma; dashed lines: partially ionized 
hydrogen plasma within a "hybrid" treatment (a model that uses the basic opacities k a from the model of partially ionized plasma and 
polarization vectors e? from the model of full ionization); light dot-dashed lines: model of full ionization. The lower curve of each type is 
related to extraordinary and the upper one to ordinary mode (the dashed curve often coincides with the solid one). 




calculated as in lPotekhin fc Chabrierl (j2Q0.1 l200l . The 
opacities that take into account partial ionization are 
plotted by solid lines, and those assuming full ionization 
by dot-dashed lines. The dashed lines correspond to a hy- 
brid approach, where the polarization of normal modes 
is described by the formulae for a fully ionized plasma, 
and k a (ui) take into account bound-bound and bound- 



free atomic transitions l)Potekhin fc Chabriertl2003l) . Al- 
though this approach is closer to reality than the model 
of full ionization, there are significant differences from 
the more accurate result drawn by the solid lines. In 
particular, the feature near hw — 0.5 keV is missed in 
the hybrid approximation. 
With increasing temperature, the differences between 
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Fig. 8. — Spectrum of hydrogen atmosphere with B = 10 12 G 
(field normal to the surface) and effective temperature T e ff = 
10 6 K. The solid line is for the self-consistent model of a partially 
ionized atmosphere, the long-dashed line is for the hybrid model, 
the short-dashed line is for a fully ionized atmosphere, and the 
dotted line is for a blackbody with T = T c g. 



the self-consistent and hybrid approximations go away. 
Figure [7| shows the case where the plasma parameters 
are the same as in Figure El except for T — 10 6 K. At 
this temperature, xn — 1.4%. Such small amount of 
neutral atoms is still very important for the opacities, 
but the hybrid approximation yields the result close to 
self-consistent one. 

Summarizing, we conclude that the hybrid approach 
to calculation of the mode opacities can be applicable if 
the fraction of bound states in the plasma, xr is small. 

6. SPECTRA 

Examples of application of the self-consistent opac- 
ity calculation for partially ionized hydrogen NS atmo- 
spheres are given in Figures 101 Here, the atmosphere 
parameters are the same as for the low-field models in 
IHo et alJ pOOl . The solid lines are obtained using the 
self-consistent calculation of the opacities, while the long- 
dashed lines reproduce the hybrid treatment described 
above. The fully ionized plasma model (short dashes) 
and blackbody (dots) are shown for comparison. The dif- 
ference in the spectra obtained using the self-consistent 
and hybrid approaches is partly due to the difference in 
the temperature profiles within the atmosphere. As ex- 
pected, this difference is small at relatively low magnetic 
field B = 10 12 G and effective temperature T eff = 10 6 K 
(Fig.|SJ), where the fraction of atoms is small at every op- 
tical depth in the atmosphere, but it becomes larger for 
lower temperature (T c g = 5 x 10 5 K, Fig.EJ) or stronger 
field {B = 10 13 G, Fig. Uni). In the case shown in Fig- 
ure El the temperature profile T(p) calculated with the 
new (self-consistent) opacities is higher by ~ 10%, so that 
photons come from shallower (lower atomic fraction) lay- 
ers of the atmosphere, which results in weaker spectral 
features due to the atomic transitions. A modification of 
the temperature profile is also responsible for the weaker 
proton cyclotron feature in this case. At contrast, in the 




ure, 



case shown in Figure ITUI the temperature profile is lower 
by < 10%, and photons come from deeper (higher atomic 
fraction) layers, resulting in stronger atomic features. 

For superstrong fields (B > 10 14 G), the current atmo- 
sphere models are less reliable because of the unsolved 
probl ems of mode c onversion and dense plasma effect 
(e.g.. iHo et al.11200.'^ . whose importance increases with 
increasing B. 

7. CONCLUSIONS AND OUTLOOK 

We have studied the polarizability and electromagnetic 
polarization modes in a partially ionized, strongly mag- 
netized hydrogen plasma. The full account of the cou- 
pling of the quantum mechanical structure of the atoms 
to their center-of-mass motion across the magnetic field 
is shown to be crucial for the correct evaluation of the 
polarization properties and opacities of the plasma. The 
self-consistent treatment of the polarizability and absorp- 
tion coefficients is ensured by use of the Kramers-Kronig 
relation. Such treatment proves to be important if the 
ionization fraction of the plasma is low (< 50%). For 
high degree of ionization (> 80%), the polarizability of 
a fully ionized plasma rem ains a good approximation , 
just as previously assumed ijPotekhin fc Chabried l2003). 
This approxima tion was adopted in th e NS atmosphere 
models built in IHo et all (|2003L |2004[) . A comparison 
with updated spectra based on the self-consistent treat- 
ment (<® shows that this approximation is satisfactory if 
B < l(P G and T cff > 10 6 K. The self-consistent treat- 
ment is needed in the atmospheres of cool or ultramag- 
netized NSs, with relatively low degrees of ionization. 

There are several limitations of the present model, 
which may become important for the magnetar fields 
and/or for low T c ff. While H atoms are treated ac- 
curately in our calculations of the EOS and opac- 
ities, H2 molecules are included in the EOS using 
the static approximation (i.e., without their CM cou- 
pling) and neglected in the opacities. Other bound 
speci es, such as H 2 + (e.g.. iTurbiner fc Lopez Vievral 
12003^ . H 3 2+ l|L6pez Vievra fc Turbinerll2000|) . and PL 
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Fig. 10. — Same as in Fig. I&lbut for higher field strength, B 
10 13 G. 



chains ijLai. Salpeter. fc Shapirolli992() are not included. 
Moreover, the NS may have a condensed surface, wit h 
negligible vapor above it ijLai fc Salpeterll9"97llLail 2001'). 
For a NS with mass M = 1.4 M© and radius R = 
10 km, estimates at logi 3 (G)= 13.5-15 (based on 
IPotekhin fc Chabrierll2004|) indicate that a thick H at- 
mosphere will be present, and a condensed surface will 
not occur, provided T cS > 3.8 x 10 5 (B/10 14 G) 1 / 4 K; 



slightly higher T ff is needed to ensure negligible abun- 
dance of molecules. 

Another uncertainty in ultramagnctized NS atmo- 
spheres is the dense plasma effect: the decoupling layer 
for photons in the atmosphere (where optical depth w 1) 
may occur at high density where the electro n plasma 
frequency exceeds the photon frequency fe.g.. IHo et alJ 
120031 lLl"ovdll2003j) . The present treatment is not appli- 
cable for such cases. In addition, construction of reliable 
atmosphere models at B > 10 14 G requires solut ion of 
the problem of mode conversion ijLai fc Holl2003aj) . 

Furthermore, for fitting observed spectra one should 
construct a grid of models with different field orientations 
and a range of field strengths, and produce angle- and 
field-integrated synthetic spectra for an assumed field ge- 
ometry. Since all the discussed spectral resonances are 
B-dependent, and some of them arc f?B-depcndent, we 
expect that such integration will somewhat smooth the 
spectral features ijHo fc Laill2004jh 
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APPENDIX 

FITTING FORMULAE FOR THE VACUUM POLARIZATION COEFFICIENTS 

The studies o f the vacuum polarization by strong fields have long history; an extensive bibliography is given by 
ISchubertJ (f2000^l . In the low -energy limit, huo <C m e c 2 , the Euler-Heisenberg effective Lagrangian can be applied. The 
relevant dimensionless magnetic-field parameter is b = B/Bg, wher e Bq = m^c 3 /eh = 4.414 x 10 13 G. In the limit 
del, the vacuum polarization coefficients are given by llAdlerlll97lD 



45tt 



To, 

45tt ' 



-^b\ 
45tt ' 



(Al) 



where ct{ = e 2 /he = 1/137.03 6 is the fine-structu r e cons tant. For arbitrary B, the tensors of vacuum polarization y vac 
and /i. have been obtained bvlHevl fc Her nouisti (Il997fl in terms of special functions and bv iKohri fc Yamadal (2002) 
numerically. The result of lHevl fc^HeniquistJ l)1997j) can be reduced to the following more convenient form: 
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where X(£) is expressed through the Gamma function T(x): 
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Results of calculation according to Eqs. 
mately represented by 



X'(0 = dX(0/d£, C = 6- 1 . (A6) 
2ll- (|A6|l agree with Fig. 2 of IKohri fc Yamadal l|2002]) and can be approxi- 
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7af 2 1 + 1.2 6 

45?r 1 + 1.33 6 + 0.56 6 2 ' 1 7 

af 6 2 



3tt 3.75 + 2.7 6 5 / 4 + 6 2 ' 



(A9) 



Equations llA7(l llA9|l exactly r ecover the weak -field limits IjAljl and the lead ing t erms in the high-field (b ^> 1) 
expansions (Eqs. [2.15]-[2.17] of lHo fc Laii r2003): in the latter regime, Eqs. O and jSHJl ensure also the terms next 
to leading. The maxim um e rrors are 1.1% at b = 0.07 for equation l|A7|l . 2.3% at b — 0.4 for equation (|A8|) . and 4.2% 
at b = 0.3 for equation (55} • 



REFERENCES 



Adler, S. L. 1971, Ann. Phys. (N.Y.), 67, 599 

Armstrong, B. M., & Nicholls, R. W. 1972, Emission, Absorption 
and Transfer of Radiation in Heated Atmospheres (Oxford: 
Pergamon) 

Becker, W. & Aschenbach, B. 2002, in 270. WE-Heraeus Seminar 
on Neutron Stars, Pulsars, and Supernova Remnants, MPE 
Report 278, ed. W. Becker, H. Lesch, & J. Triimper (Garching: 
MPE), 64 

Bulik, T., & Pavlov, G. G. 1996, ApJ, 469, 373 

Ginzburg, V. L. 1970, The Propagation of Electromagnetic Waves 

in Plasmas, 2nd ed. (London: Pergamon) 
Haberl, F. 2004, Adv. Sp. Res., 33, 638 
Heyl, J. S., & Hernquist, L. 1997, J. Phys. A, 30, 6485 
Ho, W. C. G., & Lai, D. 2001, MNRAS, 327, 1081 
Ho, W. C. G., & Lai, D. 2003, MNRAS, 338, 233 
Ho, W. C. G., & Lai, D. 2004, ApJ, 607, 420 

Ho, W. C. G., Lai, D., Potekhin, A. Y., & Chabrier, G. 2003, ApJ, 
599, 1293 

Ho, W. C. G., Lai, D., Potekhin, A. Y., & Chabrier, G. 2004, Adv. 

Sp. Res., 33, 537 
Israel, G., Mereghetti, S., & Stella, L. 2002, Mem. Soc. Astron. 

Ital., 73, 465 

Kohri, K., & Yamada, S. 2002, Phys. Rev. D, 65, 043006 

Lai, D. 2001, Rev. Mod. Phys., 73, 629 

Lai, D., & Ho, W. C. G. 2002, ApJ, 566, 373 

Lai, D., & Ho, W. C. G. 2003a, ApJ, 588, 962 

Lai, D., & Ho, W. C. G. 2003b, Phys. Rev. Lett. 91, 071101 

Lai, D., & Salpeter, E. E. 1997, ApJ, 491, 270 

Lai, D., Salpeter, E. E., & Shapiro, S. L. 1992, Phys. Rev. A, 45, 
4832 

Landau, L. D., & Lifshitz, E. M. 1989, Statistical Physics, Part 1, 

3rd Edition (Oxford: Pergamon) ^^^^^^^^^^^^ 
Lloyd, D. A. 2003, MNRAS, submitted |astro-ph/0303561] 
Lopez Vieyra, J. C, & Turbiner, A. V. 2000, Phys. Rev. A, 62, 
022510 



Meszaros, P. 1992, High-Energy Radiation from Magnetized 

Neutron Stars (Chicago: University of Chicago Press) 
Ozel, F. 2001, ApJ, 563, 276 

Pavlov, G. G., & Meszaros, P. 1993, ApJ, 416, 752 

Pavlov, G. G., & Potekhin, A. Y. 1995, ApJ, 450, 883 

Pavlov, G. G., & Zavlin, V. E. 2003, in XXI Texas Symposium 

on Relativistic Astrophysics, ed. R. Bandiera, R. Maiolino, &: F. 

Mannucci (Singapore: World Scientific), 319 
Pavlov, G. G., Shibanov, Yu. A., Zavlin, V. E., & Meyer, R. D. 

1995, in NATO ASI Ser. C 450, The Lives of the Neutron Stars, 
ed. M. A. Alpar, U. Kiziloglu, & J. van Paradijs (Dordrecht: 
Kluwer), 71 

Potekhin, A. Y. 1994, J. Phys. B, 27, 1073 
Potekhin, A. Y. 1998, J. Phys. B, 31, 49 

Potekhin, A. Y., Chabrier G., & Shibanov, Yu. A. 1999, Phys. Rev. 
E, 60, 2193 

Potekhin, A. Y., & Chabrier, G. 2003, ApJ, 585, 955 
Potekhin, A. Y., & Chabrier, G. 2004, ApJ, 600, 317 
Potekhin, A. Y., & Pavlov, G. G. 1993, ApJ, 407, 330 
Potekhin, A. Y., & Pavlov, G. G. 1997, ApJ, 483, 414 
Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 

1996, Numerical Recipes in Fortran 77: The Art of Scientific 
Computing (Cambridge: Cambridge University Press) 

Schubert, C. 2000, Nucl. Phys. B, 585, 407 
Shibanov, Yu. A., & Zavlin, V. E. 1995, Astron. Lett., 21, 3 
Shibanov, Yu. A., Pavlov, G. G., Zavlin, V. E., & Ventura, J. 1992, 
A&A, 266, 313 

Turbiner, A. V., & Lopez Vieyra, J. C. 2003, Phys. Rev. A, 68, 
012504 

Zane, S., Turolla, R., & Treves, A. 2000, ApJ, 537, 387 

Zane, S., Turolla, R., Stella, L., & Treves, A. 2001, ApJ, 560, 384 



